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o : 

, Origin of turbulence in cold accretion disks, particularly in 3D, which is expected to 

1 be hydrodynamic but not magnetohydrodynamic, is a big puzzle. While the flow must 

» I ' exhibit some turbulence in support of the transfer of mass inward and angular momentum 

' outward, according to the linear perturbation theory it should always be stable. We 

demonstrate that the 3D secondary disturbance to the primarily perturbed disk which 
^— I i exhibits elliptical vortices into the system solves the problem. This result is essentially 

' applicable to the outer region of accretion disks in active galactic nuclei where the gas 

is significantly cold and neutral in charge and the magnetic Reynolds number is smaller 
,__( | than 10 4 . 

> : 

1. Introduction 

o . 

£Nl ■ Despite much effort devoted, the origin of hydrodynamic turbulence in Keplerian 

accretion disks is still poorly understood. This is essentially important for accre- 
tion disks around quiescent cataclysmic variables, proto-planetary and star-forming 
f^i . disks, and the outer regions of disks in active galactic nuclei. 1 

| A Keplerian accretion disk flow having a very low molecular viscosity must 

Q H i generate turbulence and successively diffusive viscosity to support transfer of mass 

inward. However, theoretically this flow never exhibits any unstable mode. On the 
other hand, laboratory experiments of Taylor-Couette systems, which are similar to 
Keplerian disks, seem to indicate that although the Coriolis force delays the onset 
of turbulence, the flow is ultimately unstable to turbulence. 2 However, some other 
experiments say against it. 3 We believe that not finding hydrodynamic turbulence 
is due to their choice of large aspect ratio and small Reynolds number. 

Various kinds of secondary instability, such as the elliptical instability, are widely 
discussed as a possible route to self-sustained turbulence in linearly perturbed shear 
flows (see, e.g. 4 ). These effects have been proposed as a generic mechanism for the 
breakdown of many two-dimensional high Reynolds number flows. However, such 
effects have not been discussed properly in literatures for rotating Keplerian flows. 

Therefore, we plan to show that the three-dimensional secondary perturbation 
can generate large growth in the flow time scale and presumably trigger turbulence 
in the Keplerian disk. Unlike the two-dimensional transient growth studied earlier 5 
which were shown to be killed immediately in presence of vertical structure, in 
the present case we demonstrate essentially three-dimensional growth. Possibility 
of large growth in shear flows with rotation by a three-dimensional perturbation 
opens a new window to explain hydrodynamic turbulence which has an important 
implication not only in astrophysics but also in general physics and fluid dynamics. 
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2. Perturbation equations 

We consider a small section of the Keplerian disk in the shearing box approximation 5 
with the unperturbed velocity U p = (0, — x, 0). We also assume that the incompress- 
ible flow extends from x = — 1 to +1 with no-slip boundary conditions and is un- 
bounded with periodic boundary condition along y and z. Therefore the linearized 
Navier-Stokes equations for a 2D primary perturbation such that U% — > w x (x,y, z,t), 
If? — > U p (x) + w y (x, y, z, t), and pressure p — > p + p{x, y, z, t) are given by 

dw x dp 1 dw v ^ , „. dp 1 „, , . 

with the equation of continuity dw x / dx + dw y / dy = 0, where the angular frequency 
O = 1/q and q — 3/2 for a Keplerian disk. When the Reynolds number R e is very 
large, the solution of above equations is 

k k 
w x = £-f sm(k x x + k y y), w y = ~C-jj sin(k x x + k y y), I 2 = k 2 x + k 2 , (2) 

where C is the amplitude of vorticity perturbation. Under this primary perturbation, 
the flow velocity and pressure modify to 

U = U p + w — (w x , — x + w yi 0) = A.d, p' = p + p{x,y, z,t), (3) 

where d is the position vector and A is a tensor of rank 2. 

Now we concentrate on a further small patch of the primarily perturbed flow 
such that x << l/k x , y << l/k y . Then we consider a secondary perturbation to 
this flow such that U — > U + u and p' — > p' + p with 

(u h p) = {vi{t),p{t))ey^{ik m {t)x m ), i,m= 1,2,3. (4) 

Therefore, we obtain the evolution of a linearized secondary perturbation 

vj + A* v k + 2 e mkj n m v k = -ipkj - ^-k 2 , (5) 

kj = {-^j ) k m , k v n — k A. m v n , k — k m k . (6) 

3. Solution 

We specifically concentrate on the flow having low viscosity. Therefore, the general 
solution 6 of eqn. ([5]) can be written as a linear superposition of the Floquet modes 

v i (t)=ex P (trt)/ i (0), (7) 

where <p = mt, fi(<fr) is a periodic function having time-period T = 2n/m, and a 
is the Floquet exponent which is different at different e = (k x /l) 2 . Clearly, if a is 
positive then the system is unstable and plausibly turbulent. 

When the secondary perturbation evolves much rapidly than the primary one 
and fci(0) = (0,0, 1), a has an analytic solution given by 

o = A /Ce-(2fi3-l)(2n 3 -0- (8) 
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Clearly, for a Keplerian disk a = y%e — (4 — 3£)/9. Therefore, a Keplerian flow is 
hydrodynamically unstable under a vertical secondary perturbation if £ > 1/3. For 
other perturbations a can be computed numerically described in detail elsewhere. 6 
We now plan to quantify this by computing the corresponding turbulent viscos- 
ity. For the isotropic disk fluid, the turbulent viscosity v t = ac s h, 7 where c s and h 
are local sound speed and disk thickness respectively. On the other hand, shearing 
stress T xy =< uxu 2 >= —u t Q^- Therefore, we obtain the Shakura-Sunyaev viscosity 
parameter at a disk radius r 



For a disk with Q = 0.35, k x — 3, k y — 0.7, a vertical secondary perturbation 
evolving for time t m = 10 at a disk radius r = 15 and thickness h{r)/r = 0.05, 
a ~ 0.02 (the detailed discussion is reported elsewhere 8 ). This is interesting as a 
due to MRI computed by previous authors 9 is similar order of magnitude. 

4. Discussion 

Above results verify that the three-dimensional growth rate due to the secondary 
perturbation in a Keplerian disk could be real and positive and corresponding 
growth may be exponential and significant enough to trigger elliptical instability. 
This eventually may trigger non-linearity and then plausible turbulence in the flow 
time scale. As this growth is the result of a three-dimensional perturbation, underly- 
ing perturbation effect should survive even in presence of viscosity. We also see that 
the corresponding viscosity a to transport matter inward and angular momentum 
outward is significant and comparable to that due to MRI. 

References 

1. C. Gammie, & K. Menou, Astrophys. J. 492, L75 (1998); O. Blaes, & S. Balbus, 
Astrophys. J. 421, 163 (1994); K. Menou, & E. Quataert, Astrophys. J. 552, 204 
(2001). 

2. D. Richard, & J.-P. Zahn, Astron. Astrophys. 347, 734 (1999). 

3. H. Ji, M. J. Burin, E. Schartman, & J. Goodman, Nature 444, 343 (2006). 

4. R. Pierrehumbert, Phys. Rev. Lett. 57, 2157 (1986); B. Bayly, Phys. Rev. Lett. 57, 
2160 (1986); A. Craik, J. Fluid Mech. 198, 275 (1989). 

5. B. Mukhopadhyay, N. Afshordi, & R. Narayan, Astrophys. J. 629, 383 (2005); N. 
Afshordi, B. Mukhopadhyay, & R. Narayan, Astrophys. J. 629, 373 (2005). 

6. B. Mukhopadhyay, Astrophys. J. 653, 503 (2006). 

7. N. Shakura, & R. Sunyaev, Astron. Astrophys. 24, 337 (1973). 

8. K. Saha, & B. Mukhopadhyay, "Turbulent viscosity in the Keplerian accretion disks 
due to hydrodynamic instability" (submitted). 

9. A. Brandenburg, A. Nordlund, R. F. Stein, & U. Torkelsson, Astrophys. J. 458, L45 
(1996); J. M. Stone, J. F. Hawley, C. F. Gammie, & S. A. Balbus, Astrophys. J. 463, 
656 (1996). J. C. B. Papaloizou, & R. P. Nelson, Mon. Not. Roy. Astron. Soc. 339, 
983 (2003). 



a = — 



qn 2 (^) 3 Mr 2 ' 




(9) 



